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ABSTRACT 

Cosmological simulations of the low-density intergalactic medium exhibit a strikingly 
tight power-law relation between temperature and density that holds over two decades 
in density. It is found that this relation should roughly apply Az ^1 — 2 after a 
reionization event, and this limiting behavior has motivated the power-law parame- 
terizations used in most analyses of the Lya forest. This relation has been explained 
by using equations linearized in the baryonic overdensity (which does not address why 
a tight power-law relation holds over two decades in density) or by equating the pho¬ 
toheating rate with the cooling rate from cosmological expansion (which we show is 
incorrect). Previous explanations also did not address why recombination cooling and 
Compton cooling off of the cosmic microwave background, which are never negligible, 
do not alter the character of this relation. We provide an understanding for why a 
tight power-law relation arises for unshocked gas at all densities for which collisional 
cooling is unimportant. We also use our results to comment on (1) how quickly fluctu¬ 
ations in temperature redshift away after reionization processes, (2) how much shock 
heating occurs in the low-density intergalactic medium, and (3) how the temperatures 
of collapsing gas parcels evolve. 

Key words: cosmology: theory - cosmology: large-scale structure - quasars: absorp¬ 
tion lines - intergalactic medium 


1 INTRODUCTION 


Cosmological simulations of the low-density intergalactic 
medium exhibit a tight power-law relation between tem¬ 
perature and density that hol ds over two decades in den¬ 
sity (e.g. iHui fc Gnedinl Il997l : see Fig. [I|. This relation 
is established Az ~ 1 — 2 after the end of cosmologi¬ 
cal r e ionization processes (Hui fc Gnedinl 1 19971 : Trac et ^ 
I2OO8I : iFurlanetto &: OhI l2009l :~ McQuinn et al.l 20091 ) . Moti¬ 
vated by this limiting behavior, a power-law parameteriza¬ 
tion for the thermal state is adopted in almost all analyses of 
the Lyg forest ( e.g., [McDonald et ahlboosl : IViel et al]l2006l : 
Ibidz et aDl201Cll )n Here we provide a simple explanation for 
the physics that sets the asymptotic relation as well as an 
understanding for how quickly this asymptote is reached. 

While the mathematics of what sets the asymptotic 
temperature and its power- law slope in density w ere worked 
out in the seminal study of iHui fc Gnedinl (1 19971 ). their lin¬ 
earized approach does not explain why the same power- 
law holds over a range of ~ 100 in density. Subsequently, 
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^ Most studies that use numerical simulations implicitly adopt 
such a parametrization, as standard cosmological simulations em¬ 
ploy uniform radiation backgrounds - an approximation which 
guarantees a near power-law relation. 


iTheuns et al.l ill 9981 ) and iPuchwein et al.l ll2014l) attempted 
to provide more intuitive derivations that we show are prob¬ 
lematic. In addition, previous studies did not address how 
non-adiabatic cooling processes affect this relation, namely 
Gompton cooling off of the cosmic microwave background 
(GMB) and recombination cooling. Both cooling processes 
can be many tens of percent of the photoheating rate at 
the mean density of the Universe (and comparable to pho- 
toheating at other relevant densities). 

We show that the power-law form of the T — A rela¬ 
tion occurs because the photoheating rate per happens 
to be nearly independent of nt,, where ni is the number 
density of baryons, driving all cosmic gas towards a single 
adiabat. This form for the photoheating rate also results in 
the temperature of a gas parcel being essentially indepen¬ 
dent of its collapse history - explaining why the asymptotic 
T — A relation holds so tightly despite the varied histories 
of cosmological fluid elements. (The tightness of this rela¬ 
tion is absolutely striking in calculations that do not allow 
for shocking but allow for varied collapse histories. See the 
Zeldovich approximation calculations in Fig. m) We further 
show that Gompton cooling off of the GMB has no effect 
on the asymptotic slope of T(A) and that recombination 
cooling is just weak enough not to break the relation. 

This paper is organized as follows. We first enumerate 
the relevant heating and cooling processes in Next, in ij3l 
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2 M. McQuinn & M. Zaldarriaga 


we derive an analytic expression for the asymptotic T — A 
relation in the limit of only photoheating and adiabatic heat¬ 
ing/cooling. These are the processes considered in previous 
attempts to nnderstand this relation analytically. Section [4] 
generalizes our solutions to include all relevant cooling pro¬ 
cesses. We briefly consider three applications of our results 
in !j5l ( 1 ) how quickly temperature fluctuations are erased, 
(2) the importance of shock heating, and (3) the thermal 
evolution of collapsing gas clouds. Lastly, JHl contrasts our 
results with previous attempts to understand the intergalac- 
tic T — A relation. 

Throughout we use the standard phenomenological 
parametrization 

T = ToA^-^ (1) 

for the temperature of the low-density IGM, where A is 
the baryonic density in units of its cosmic mean. To is the 
gas kinetic temperature at A = 1, and 7 — 1 is a power- 
law index. (Sometimes 7 is called the “equation of state.”) 
When necessary, we assume a fiat ACDM Universe with 
h = 0.7, flm = 0.3, Q b = 0.045, and Yue = 0.25. Following 
iHui fc GnedinI (Il997ll . we use the Zeldovich approximation 
to follow the density (necessary for computing the tempera¬ 
ture) of fluid elements for certain calculations. Our Zeldovich 
approximation calculations draw matter elements randomly 
from the probability distribution of collapse eigenvalues, as¬ 
suming that the gas traces the dark matter. This probabil¬ 
ity distribution only depends on the variance of the density 
field, which we set to 1.2. This choice is equal to the variance 
of a real space top hat filter containing a Lagrangian mass 
of 10^° Mq in our reference cosmology at 2 ; = 3 (and ap¬ 
proximately the Jeans’ mass near the cosmic mean density 
at this redshift)0 In addition, we also use a 20 Mpc, 512® 
gas and 512® dark matter particle simulation that was run 
with the GAD GET-3 smooth particle hydrodynamic code 
llSpringelll2005l L Reionization occurs instantaneously in this 
simulation at 2 : = 9, heating all gas to 20, 000 K, and the 
subsequent heating and cooling takes rates applicable for 
primordial gas if only the H l and He l are kept ionized with 
an ultraviolet background that is flat in specific intensity 
[ergs“® Hz“® cm“^ sr“^]. 


2 RELEVANT HEATING AND COOLING 
PROCESSES 

The evolution of the temperature of a Lagrangian fluid ele¬ 
ment in an FRLW background cosmology follows 

^ = _2HT +(2) 
dt 3A dt SksTib dt ’ 

where A is the baryonic density in units of the cosmic mean 
and rib is the number density of all “baryonic” particles in 
the plasma (including electrons). Equation Q is valid when 

^ We do not omit regions that have gone through a singularity. 
Regions in which an odd number of dimensions have collapsed 
result in negative densities that are not seen in our plots. If we 
were solving the full nonlinear theory, parcels that collapsed along 
one or more dimensions typically would result in A > 10, roughly 
above the densities considered here. 



Figure 1. Temperature-density relation in a cosmological simu¬ 
lation (showing 300 randomly selected locations) and in the Zel¬ 
dovich approximation calculations (following 300 Lagrangian el¬ 
ements). Both the simulation and the Zeldovich approximation 
calculations are initialized with To = 20, 000 K and 7 — 1 = 0 at 
z = 9, and are evolved forward with Obk = 0. The resulting (T, A) 
values fall on almost a single power-law relation over two decades 
in A, especially in the Ze ldovich approx i matio n calculation. A 
similar result was found in iHui Sz Gnedir] lll997lV The simulation 
points show more dispersion than the Zeldovich approximation 
ones because of (mostly weak) shocking f ^5.2ll . 



Figure 2. Importance of different processes in shaping the tem¬ 
perature at the cosmic mean density, Tq (top panel), and the 
power-law index of the T — A relation, 7 — 1 (bottom panel). 
These calculations are initialized with (To, 7 — 1) = (20, 000 K, 0) 
at 2 : = 9. They assume that Obk = 0 to calculate the photoheat¬ 
ing rates and that the He ii has not been reionized such that the 
gas is comprised primarily of H ii and He ii, with trace amounts 
of neutrals. To evaluate 7 — 1, we assume that A = 1 and that 
growth scales so A — 1 oc a. 
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On the intergalactic temperature-density relation 3 


Table 1. Heating and cooling rates per baryon for a photoionized IGM, in units of 3860 K per free particle per Gyr.* The amplitude 
and power-law indices are evaluated at r 4 = 1, A = 1, and Z 3 = 1, where Z 3 = (I + z)/4 and T 4 = T/10^ K. The power-law indices in 
Z 3 and A are exact, and the power law in T 4 is accurate over relevant temperatures. 


H I photo. 

He I photo. 

He II photo. 

H II recomb. 

He III recomb. 

Compton 

Free-free 

Hubble 


rp—Q.7 ^3 A 
n 1 0 -'4 ^ 3 ^ 


- 0.11 Tf'2z| A 

-0.20Tf'3z|A 

-0.28 r4Z| 

-0.05 yTTZ|A 

-1.6T4Z|/^ 

l + “bk/2 

l + /{cbk) 

l+«bk/2 


‘All rates assume that the applicable bound states have ionization fractions that are either nearly zero or nearly one, showing the 
choice that maximizes the rate. We have omitted collisional cooling and He II recombination cooling, the former is subdominant at 
A < 10 and the latter is always small, with a rate of —0.011T|’ ®Z| A. The photoheating rates assume photoionization equilibrium, and 
the normalization of the H I and He I rates are calculated for a background with = 0 until an abrupt cutoff at 4 Ry. See footnote [3] 

for additional details such as the definition of /(obk)- 



Figure 3. Same as Figure [2] except for two differences: First, 
the calculations are initialized with (To, 7 — 1) = (15,000 K,0) 
at 5 : = 3.5. Second, we take photoheating and cooling rates that 
assume the He II has been reionized. 

the number of particles remains fixed - an excellent approxi¬ 
mation after reionization processes. Table[T]lists the different 
primordial gas heating and cooling processes that contribute 
additively to the {3kBni,/2)~^dQ/dt term in equation Q, 
aside from the column labeled “Hubble” which corresponds 
to the —2HT term. All of the listed processes are well ap¬ 
proximated as power laws in = T/10"^K, and the power- 
laws in A and Z 3 = (1 -|- 2)/4 are exact. The only cooling 
processes that are not included in Table [T] are those due to 
collisional excitations, which become important at A > 10 
(see 115.311 . 

Photoheating is the only non-adiabatic heating source 
included in our calculations and in standard models for 
the thermal history. These models appear to be sufficient 
to describe intergalactic temperature observations over the 

measured range of 1.5 < _ z < 5 dPuchwein et al.l l2014l : 

lUpton Sanderbeck et al.ll2015l l. Photoheating inputs an en¬ 
ergy per H I ionization of eni ~ i?ion,Hi/(2 -|- Obk) before 
He II reionization (dominated by H l photoheating; Sion ,hi = 
13.6 eV), and eni ~ 3 Eion,Hi/(2 -I-Obk) after (dominated by 
H I and He ll photoheating; see Table[T])13 Nicely, the photo- 

® Our H I and He I photoheati ng rates assume an abru pt cutoff 
at 4 Ry as found in models (e.g. lHaardt fc Madauir2012l ~). For the 


heating rates, which for a highly photoionized IGM are ex 
times the recombination rate, only depend on the spectral 
index of the background specific intensity near the ionization 
edge of the species in question (and not on the amplitude of 
the background). We parametrize the background specific 
intensity as Jv oc [ergs“^ Hz~^ cm“^ sr“^]. Ionizing 

background models predict Obk « 0 — 1 over « 1 — 10 Ry, 
with this value being determined by the intrinsic spectrum 
of the sources and the reprocessing fr om intergalactic ab¬ 
sorptions (e.g.. IHaardt fc Madaul [20121 1. This range in Obk 
translates to a 30% uncertainty in the post-reionization pho¬ 
toheating rates. We choose Qbk = 0 for the calculations in 
this paper. 

Figure [ 2 ] illustrates how Tq and 7 — 1 evolve. Here, 7 — 1 
is evaluated at A = 1; we will show that a single power-law 
holds at all densities only in the limit of late times. These 
calculations are initialized with Tq = 20, 000 K and 7—1 = 0 
at 2: = 9, consistent with simul ations and analytic cal¬ 
culat i ons of hydrogen reionization llMiralda-Escude fc ReesI 
1 1994 iTrac et 'SII 2 OO 8 I : lMcQuinnll2012l l. and use the heating 
and cooling rates that apply before He ll was reionized. The 
solid curve includes all processes that affect To, whereas the 
other curves turn off the specified process (es) to illustrate 
how they impact the thermal history. The main source of 
cooling at A ~ 1 is the expansion of the Universe. Comp¬ 
ton cooling and recombination cooling are smaller but never 
negligible, with Compton cooling being quite important for 
shaping Tq. 

Figure |3] is the same as Eigure [2] except that it considers 
the case starting with Tq = 15,000 K and 7 — 1 = 0 at 
2 = 3.5, v alues more consistent with the IGM after H e ll was 
reionized (|McQuinn et al.ll2009l : iBecker et al.l 1201 ill , and it 
takes the applicable rates for after He ll reionization. In this 
case, Compton cooling is less important and neglecting non- 
adiabatic cooling processes is a much better approximation. 


He I photoheating rate listed in Table [T] /(obk) some function 
of “bk which equals zero for Qbk = 0 and would equal Qbk if 
there were no cut-off in the spectrum at 4 Ry. For H I, there 
should be a similar functional form that depends on the cutoff at 
Tmax = 4 Ry; the Qbk scaling included in the table approximates 
this function with its Fmax —t 00 form (which is accurate to 10 %). 
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4 M. McQuinn & M. Zaldarriaga 


3 DERIVATION OF T - A RELATION 

IGNORING NON-ADIABATIC COOLING 


and parametrizing the density evolution of a gas parcel as 
A(a') = A(a) x (a'/a)", equation ([9ll becomes 


Our aim is an analytic understanding of why the temper¬ 
ature is driven to a single T — A relation. Let us intro¬ 
duce the adiabatic variable rj = T/rdj^^, using that uh ~ 
h-rf.o A(a), where nu is the total number of hydrogen 
(H I -b H ll), subscript “0” indicates the present day value, 
and A(o) is the density contrast with respect to the mean. 
The time derivative of r) is 


2/3^77 dT 2T dA 


(3) 


such that, after substituting this into equation (H) , our equa¬ 
tion for dT/dt becomes 


2/3 d^ _ 2 dQ 

^ dt SksUb dt 


(4) 


We first specialize to the simplified case of an IGM with 
only phot oheating and a d iabat ic processes, the limit con¬ 
sidered in iHui fc GnedirJ (Il99i1 i. While we saw in 321 that 
this limit is not a good approximation for the evolution of 
To at high redshifts owing to Compton cooling off of the 
CMB, this limit is a better approximation for the evolution 
of 7 — 1 (Fig. [^. In addition, for calculations that start at 
lower redshifts, this photoheating-only limit is an excellent 
approximation for both Tq and 7 — 1 (Fig.[3|. With only the 
photoheating contribution to dQ/dt, equation (|4]) becomes 


2/3 dry 2eHirHinHi 

^ = Sksub ’ 

where Fhi is the H l photoionization rate. Assuming that 
the gas is in photoionization equilibrium so that nni ~ 
QAUenif/rHi and changing the time variable to the scale 
factor, a, yields 

2 / 3^77 eHiOAUira^^^ 


Note that to reach this equation we assumed Tie ~ tijj 
and H — , appropriate in the assumed cos¬ 

mology at 2 : > 1. In what follows, we use that the CASE 
A recombination coefhcient can be written as aA(T) = 
QA(rp) X {T/Tp)-°-'^ = aA{Tp) X {n]i^n/Tp)-°-\ where 
Tp is some reference temperature. Since it makes the al¬ 
gebra cleaner, we will also approximate —0.7 as —2/3, but 
will comment on the generalization. Integrating equation 
over time yields 

rd77 77=/^ = A (7) 

J r}i o-i 

A = aA{Tp)T^2^ui ^ 

skbHonU 

or 

,5/3 _ ^5/3 ^ £ da'a'^/^^Aia')-^/^. (9) 


Since the density of most intergalactic gas decreases con¬ 
siderably over a cosmologically significant period, whereas 
the temperature decreases only by as much as a factor 
of ~ 2 , ( 777 / 77 )®^® decreases strongly with increasing time, 
erasing memory of the initial conditions as ~ ( 07 / 0 )^°^® 
for gas at the mean density. Taking the limit 77 777 


r(a)®/® 


5A7ig,oa~®/^A(fl) / r^iii/6-"/9 

3 (11/6 — 7i/9) \ I. a J 


.( 10 ) 


Equations and uni) demonstrate that the asymp¬ 
totic temperature of a gas element depends very weakly on 
whether it collapses to its current density, A (a), quickly or 
slowly. The time dependence of A(a') only enters in the tem¬ 
poral integral on the RHS of equation and then to the 
— 1/9 power. Equation (1101) shows that unless |n| > 9, the 
collapse/expansion has less than a factor of two effect on the 
proportionality constant in the asymptotic T — A relation¬ 
ship relative to the n — 0 case. The exponent n = 9 corre¬ 
sponds to A(t) changing by a factor of two over a timescale 
of da/a = 0.1. Since Aa/a = HAt, and the dynamical time 
is for A 1, Aa/a ss 0.1 corresponds to the 

dynamical time for gas with A « 100. Using these relations, 
the expected n for A = 10 gas (the highest densities of inter¬ 
est) is « 3, which would lead to a 20% deviation from the 
71 = 0 normalization. Within trajectories in the Zeldovich 
approximation (where the collapse happens too suddenly at 
high densities), we find visually a range of 5 < n < 7 for 
A = 10 and 2 = 30 with n = 5 resulting in 15% smaller 
temperatures than n = 7 in equation cni). In detail, the 
Zeldovich calculations shown in the bottom panel of Fig¬ 
ure [2] have a standard deviation in AT/T of 5% at A = 10, 
and this becomes somewhat smaller at lower A. Thus, the 
tightness of the T — A relation owes largely to the weak de¬ 
pendence of the final temperature on the previous collapse or 
expansion history of a gas element. The average dependence 
of 71 on A does induce a small excess slope to the T — A 
relationship that we estimate in 34] 

In the 71 = 0 case, equation m becomes 


T{a,A) 

To,iim(a) 


To,iMA^/^ 

/ 10aA(Tp)T^''^eHinH,oa-^^^(l - [^]^^^*^) 

V sskiHonU^ 


/ EHI 

3/5 


V5eV/ 


0.25 

'v J 


although we emphasize that the result changes little for any 
71 that is likely to apply in the low density IGM. Note again 
that CHI ~ 13.6eV/(2-b Obk) before He ll was reionized, and 
it equals three times this value afterwards. 

In conclusion, in the photoheating-only limit we Hnd 
T oc A3/5, 

in excellent agreement with the low-redshift 
asymptote of the solution to T{z,A) with all cooling pro¬ 
cesses included. (If we had instead used the parametriza- 
tion oa ~ T~^ for the recombination rate, we would find 
T oc with 7 — 1 = 0.59 for /i = 0.7.) The physical 

reason for this asymptotic relation is quite apparent from 
the form of equation ®: 


Qphoto(^! ) 


( 12 ) 


This is the effective n that corresponds to the instantaneous 
evolution of dA/dt. 
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where (/photo(??, riir) is the heating rate of the gas. Because 
%hoto/n^if^ is almost independent on uh at fixed 77 at 2 < 6 
(scaling as 77 “°'^n 5 f^), all gas is driven towards the same 
adiabat, erasing memory of the initial conditions. Because 
the heating is nearly constant along an adiabat, this con¬ 
stancy also implies that the final temperature of a gas parcel 
depends extremely weakly on its collapse history. We now 
show that these conclusions hold when all relevant heating 
and cooling processes are included. 


4 THE RELATION WITH ALL COOLING 
PROCESSES INCLUDED 


So far we have considered the case of photoheating and 
only adiabatic cooling processes. Taking r 4 = 1, A = 1, 
Z 3 = 1 , and able = 0 , the amount of atomic and free-free 
cooling is about 50% of the amount of photoheating before 
He II was reionized, and about 20% afterwards (see Table[T]). 
Larger percentiles occur for other parameter choices. Thus, 
it is plausible that cooling processes break the tightness or 
power-law form of the T — A relation at certain epochs and 
certain densities. This section first provides an initial sketch 
of why cooling processes do not break the relation. Then, in 
H.ll we treat the case of Compton cooling in more detail. 

Let us treat the non-adiabatic coolants as a small cor¬ 
rection to our previous solution such that rj = -|- 

where is the photoheating-only solution and is the 
perturbation around this solution including all other pro¬ 
cesses. We also multiply all cooling terms that are not in 
our “zeroth order” equation, equation ©, by e. Collecting 
the terms that are first order in e yields the equation 


2 / 3 ^ (0) 

H ^ 4 ^ ' yphoto ^(0) 


dt 


(^Smp + , (13) 


where we use the notation that is the heating/cooling 
rate of X [degrees Kelvin per second per particle] evaluated 
at 77 *-°^ and tih, and where we have grouped all recombi¬ 
nation cooling rates into qrec- To derive a rough estimate 
for the size of the correction, we further assume that the 
term is smaller than the others so that this equa¬ 
tion becomes 


0.7q 


(0) ^ 
photo ^(0) 




^Comp 


+ 




(14) 


We find numerically that this equilibrium solution is a de¬ 
cent approximation for the full solution to equation m- 
Using Table[T]to insert physical values for the heating rates, 
equation m becomes 


y(i) 

r( 0 ) 


recombination+ff 


Compton 

- a3^-0.15A°“ 



(15) 


before the He ll was reionized and for a = 0. To represent 
the formula in this simple form, we assumed that the recom¬ 
bination cooling and free-free cooling rates scale as bi¬ 
secting the scaling of these different rates. If we had instead 
substituted the applicable rates for after He ii reionization, 
the fractional contribution of recombination plus free-free 
cooling would stay nearly the same, while the contribution 
of Compton cooling would be reduced by a factor of nearly 
three relative to equation m- 


Equation GSl) shows that Compton cooling does not al¬ 
ter the T — A relation. This occurs simply because 

with no direct A dependence. Recombination cooling 
and free-free cooling do alter the T— A, resulting in its frac¬ 
tional contribution to having a strong dependence on 
A, as A'^ ® in equation (O. However, even at the highest 
densities where the T — A relation applies, A ~ 10, recom¬ 
bination and free-free cooling contribute a significant 40% 
correction for = 7000 K. Using Equation (I15II . the 

correction from recombination and Compton cooling to the 
equation of state results in a flatter value of 


7-1 


recombination+ff 


avg. collapse history 


0.59- 


0.06 


1 - 0.35^3 


T 


(0) 


0 


7000K 


-b 


0-03 ^ 0.5 

1 - O. 35 Z 3 


(16) 

assuming Qbk = 0 and using the P = 0.7 value for the lead¬ 
ing order rather than our simplification of )? = 2/3. (Ignore 
the “avg. collapse history” term in this expression momen¬ 
tarily.) A similar correction to the “recombination-!-ff” term 
in the above equations is also noted seen in the full calcula¬ 
tion of T{z)-. Compare the low-redshift asymptote of the “all 
processes” curve to the “no recombination cooling curve” in 
the bottom panel of Figure [21 

The “avg. collapse history” term in equation (O is the 
correction from adiabatic heating (the n-dependent term in 
To™). In particular, we assumed that this correction scales 
as A®’® as expected if the characteristic time is the dynam¬ 
ical time, with n = 3 corresponding to A = 10 and n = Q 
to A = 1 as argued for in llsli This correction acts to cancel 
half of the “recombination+ff” term. 

The above assumed that Compton cooling can be 
treated perturbatively, which Figure [2] shows is a dubious 
assumption, especially at high redshifts (although it looks 
to be a good assumption for the case shown in Fig. 121 ). We 
show in M.ll that Compton cooling also does not affect the 
asymptotic 7 — 1 even when this process shapes the thermal 
history. 


4.1 Full treatment of Compton cooling 

Here we calculate r(A, a) in the presence of Compton cool¬ 
ing. Unlike the previous section, we do not make the dubious 
assumption that Compton cooling contributes a small cor¬ 
rection to the full solution. The evolution of the temperature 
for the case of photoheating plus Compton cooling can be 
solved exactly. The differential equation for the temperature 
becomes (after a few simple manipulations) 


3^ 

5 da 




Aad^'^n 


-(2/3)^ + l-2/3 
H I 


where 


4(JT Ur ad , 0 

SrriecHoQm^ 


0 . 011 , 


(17) 


(18) 


^ We note that the “avg. collapse history” correction is exactly 
what is needed to explain the 7 = 0.6 2 asymptote derived using 
linear theory in iHui & Gnedii] l|l997l'l for the photoheating-only 
limit. 
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6 M. McQuinn & M. Zaldarriaga 


and flrad is the radiation constant. We solve 


^^+Ca-’^^^G = 5{a-a), (19) 

for the Green’s function G of the linear in 77 ®^^ part of equa¬ 
tion GZl, which yields 


G{a, a) 



^9{a — a), 


( 20 ) 


where 9 is the Heaviside function. The solution to equa¬ 
tion (II3 is the sum of a term that satisfies the linear in 
part of this equation plus the integral of G over the 
source term: 




,-V2_„-5/2 


Q® f da a 
J a,,- 


5/6 _2c(a'-5/2-a- = / 2 ) 


A{a'y 


( 21 ) 

-1/9 


Notice that our equation is identical to equation © in the 
limit of no Compton cooling (G —^ 0). Our s olution also re¬ 
duces to the solution for To in Appendix A of iMallov &: Lidj 
(|2niRh for the case A(a) = 1. At high redshifts, Compton 
cooling enhances the rate that memory of the initial con- 
ditions is erased via the factor e~ ““ >. Equa¬ 

tion m also shows that once the memory of the initial con¬ 
ditions is erased, 17 ®/® oc or equivalently T oc A^^®, 

the same asymptote that we found in the photoheating-only 
limit. Also as in the photoheating-only limit, the history of 
the density contrast enters in the temporal integral to the 
— 1/9 power, again making the result insensitive to the tem¬ 
poral evolution of A(a). 

We find that Compton cooling results in an extra factor 
in the asymptotic relation for T{A) over our photoheating- 
only solution (eqn. Hill : 


^ 15/lle-(^)"^'a-“/®(aii/6_a,ii/6) ^ 


where 

-^Comp 


r 11 


11 


L 15- 

V a / J 

15 ’ 

1 J J 


T(a,A) = 


Comp 


T' 
Tb.lim 


where T is the incomplete gamma function, flcrit = 
[2G/3]^^® = 0.14, and we have assumed n = 0. Evaluat¬ 
ing Acomp numerically yields {0.81, 0.78, 0.76} at 2 = 
(3,4, 5} for Zi = 9. For Zi = 6 , Acomp instead yields 
{0.88, 0.89, 0.92}. These small reductions in the asymptotic 
To are not the reason Compton cooling has such a large ef¬ 
fect in Figure 2; the primary effect of Compton cooling is in 
its exponential suppression for a > Ocrit of Ti. 

In conclusion, while the asymptotic To is marginally 
altered by Compton cooling off of the CMB, the asymptotic 
7 — 1 is unaffected by it. Compton cooling does not increase 
the sensitivity to a gas parcel’s collapse/expansion history, 
and it accelerates the convergence to the asymptotic T{A). 


5 APPLICATIONS OF RESULTS 

Here we use our results to understand (1) how quickly tem¬ 
perature fluctuations are erased 1 45.111 . (2) the importance 
of shock heating (SEU, and (3) the thermal evolution of 
collapsing gas clouds 1 45.311 . 


5.1 timescale to reach the asymptote and erase 
temperature fluctuations 


Let us consider a toy universe that has an isothermal T — A 
relation after reionization. The temperature at a later time 
follows from equation HSU: 


j.5/3 


rrpComp-iS/S 

1-^0,11111 J 


A-f 


a~^A 


C = e 


5/2/ —5/2 —5/2\ 


10/9 


( 22 ) 

(23) 


where Zcrit = l/flcrit — 1 = 6.15 and Ti is the temperature 
after reionization. The timescale for the T—dependent term 
to become a faction / of is 

a/a, (24) 

where we have approximated the 10/9 exponent as unity. 
For example, if Ti is twice - roughly the multi¬ 

ple expected after H i and He ii reionization -, it takes 
a/tti « 1.9(C/Ai)^/® for the temperature to equal half of 
its limiting value (/ = 0.5). Interestingly, this timescale is 
slightly shorter in initially overdense regions and there is no 
dependence on the final density. 

Fluctuations in temperature - an inevitable byproduct 
of a temporally extended and inhomogeneous reionization 
process - redshift away in a nearly identical manner to devi¬ 
ations from the asymptotic temperature. After reionization, 
the spatial variance of 77 ®/® is simply 

var[77®/®]A = C^var[77®^®]A, (25) 


where var[y]A is the variance of Y fixing the final density 
A but averaging over all spatial positions, and the subscript 
i again denotes some initial time. The fractional variance in 
the temperature in terms of the mean is then 


T 


/i2 

' C var 


Vi 


iv) 


(26) 


specializing to times when ATjT <C 1 so that var[y"*] « 
?Ti^var[y]. This approximation likely holds soon after reion¬ 
ization processes as simulations show th at AT jT ^ 1 at 
the end of H l and He ll reionization jTrac et al.l l2008l : 
iMcQuinn et al.l I2OO9I : ICompostella et ^ l2013l i. Thus, the 
RMS of temperature fluctuations decreases over essentially 
the same timescale the average temperature asymptotes, 
again with a slightly quicker relaxation for initially denser 
regions. 

The above discussion relates the temperature fluctu¬ 
ations generated at some “initial” time to those at some 
time thereafter. A better way to conceptualize how tem¬ 
perature fluctuations are imprinted, particularly during hy¬ 
drogen reionization, is that a region is heated to a simi¬ 
lar temperature when it is reionized but different regions 
are reionized at different times. The duration of reionization 
and the associated cooling results in a spre ad in temperature 
l|Trac et al.l I2OO8I : iFurlanetto h OhI l2009li . In this scenario, 
the above relations can be used to calculate how the tem¬ 
perature fluctuations relate to the duration of reionization. 


5.2 the importance of shocks 

We have shown that without shocks a very tight T — A rela¬ 
tion should hold at A < 10 (see Fig. 1). Thus, in a cosmolog¬ 
ical simulation, any deviations of the fluid elements from a 
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Figure 4. Fraction of gas elements in our cosmological simulation 
that are shock heated by a factor of 1.5, 2 and 10 over the mean 
T(A) (thick, medium width, and thin curves, respectively) at 
2 = 2.3 (solid curves) and 2 = 3.0 (dashed curves). The analytic 
calculations in this paper are valid to the extent shock heating is 
not important. 


tight T — A relation indicates shocking. The black points in 
Figure 1 show the T{A) of 300 particles randomly selected 
from a cosmological simulation. The particles do not show 
as tight a T — A relation as in the Zeldovich approximation 
calculations (red points), suggesting that some of the gas 
elements in the simulation have experienced weak shocking, 
although most particles fall within 10% of the expected re¬ 
lation. Figure [4] shows the fraction of gas that is heated by 
a factor of 1.5, 2 and 10 over the mean T{A) as a function 
of A (thick, medium width, and thin curves, respectively) 
at 2 = 2.3 and 2 = 3. By 2 = 2.3, the simulation shows 
that only a few percent of fluid elements with A = 1 have 
been shock heated significantly, with this fraction scaling as 
~ A° ®. The shocked fraction in the simulation decreases 
quickly with increasing redshift. 


5.3 the thermal history of fluid elements 

The maximum Jeans’ mass of a collapsing matter cloud, and 
whether this mass is larger than the cloud’s gravitational 
mass, determines whethe r it can overcome pressu re and col¬ 
lapse on to a galaxy (e.g. iNoh fc McQuinnIuOldl) . Since the 
Jeans’ mass depends on the gas temperature, the thermal 
history of collapsing fluid elements sets the minimum mass 
of galaxies. The solid curves in Figure [5] show the thermal 
history of unenriched gas that collapses at the specified red- 
shifts. These curves are calculated assuming spherical top 
hat collapse, that the gas was heated to 20, 000 K at 2 = 9 
by the reionization of H l and He l, that it was kept ionized 
with a background of Obk ~ 0 thereafter, and that there 
was an instantaneous reionization of He ll at 2 = 3 (and an 
associated heat injection of 7000 K). While modeling reion¬ 
ization processes as instantaneous is artificial, the thermal 
history of this scenario roughly approximates the measured 



Figure 5. Solid curves show the T(A) trajectories of gas parcels 
that collapse at the specified redshifts, assuming spherical top hat 
collapse and assuming that the gas was heated by reionization to 
20, 000 K at 2 = 9. Gas parcels start at this temperature and 
evolve with time in a clockwise fashion in this plot. The dashed 
curves show the T — A relation at the specified redshifts, with 
the filled dot on each curve indicating the temperature at the 
cosmic mean density. These curves show that a power-law relation 
generally holds about a decade in density above the cosmic mean, 
after which collisional cooling processes become important. These 
calculations also assume a heat injection of 7000 K at z = 3 to 
approximate the heating from He II reionization. 


history (iBecker et akllioilh . The trajectory of a fluid ele¬ 
ment is clockwise in this plot, with uh c>o as z ^ -^^iiiH 
The dashed curves in Figure [5] show the T — A relation 
at the specified redshifts, with the filled dot on each curve in¬ 
dicating the temperature at the cosmic mean density. These 
curves can be used as a clock to infer the redshifts of dif¬ 
ferent points along the (solid) collapsing curves. They also 
show that a power-law T — A relation generally holds about 
a decade in density above the cosmic mean, to densities at 
which collisional cooling processes become importantlj| The 
solid curves in Figure [S] resemble the qualitative picture for 
a collapsi ng fluid element’s trajec tory in T — A space pre¬ 
sented in iNoh fc McQuinrJ ll20l4 l. which this study then 
used to infer the minimum mass of galaxies. 


® Of minor importance to our results, we have also assumed that 
at densities where collisional cooling processes become important, 
th e He ll is self-shie l ding to ionizing backgrounds, as motivated 
in iNoh &: McOuinnI ll20l4l . so that He II collisional cooling can 
occur. 

^ Perhaps surprisingly, the asymptote of the solid curves bifur¬ 
cates into two values as njj —>■ oo. We believe this owes to the 
two peaks in the cooling function of primordial gas as a function 
of temperature (owing respectively to H I and He II cooling), 
resulting in two thermally stable equilibria. There is evidence 
that a similar bifurcat ion happens in cosmological simulations 
dNoh &: McOuinr]|2014h . 
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6 PREVIOUS EXPLANATIONS FOR 7 

We are aware of two previous attempts to pr ovide an intu¬ 
itive e xplanation for the T —A relation. First. [Theuns et al.l 
(Il998l . Appendix C) explained the relation as owing to the 
balance between photoheating and adiabati c cool ing from 
the uniform Hubble expansion. iTheuns et al.l lll998l l derived 
the equation 7 — 1 = l/(/3 + 1) with /3 « 0.7, the same re¬ 
lation derived in SJS] However, we think this approach gives 
the correct result from an incorrect insight, as the asymp¬ 
totic 7 — 1 has nothing to with cosmological expansion as 
this approach posits. Our derivation shows that photoheated 
gas that cannot cool converges to the same 7 — 1 , virtually 
independent of its collapse or expans ion history. We offer 
more quantitative musings about the iTheuns et al.l ill998l l 
app roach in the ensuing fo otnote^ 

IPuchwein et al.l (l2014l l provided a separate explanation 
for 7 — 1. They derived 7 — 1 by an argument that is 
equivalent to taking the t = 00 solution to equation 
for the case of constant A(t). Thus, they solved dT/dt = 
(3/2n6fci,)“^ [d(5/dt]photo in the limit of large times. The so¬ 
lution to this e quation yields the corr ect slope of 7 — 1 = 
1/(1 -I- j3). The IPuchwein et al.l il2014l ~) approach does hold 
insight into what sets the asymptotic slope, and indeed is 
our solution specializing to the case in which the density of 
a gas parcel is constant in time (our solutions for n = 3). 
However, this explanation ignores the adiabatic cooling as¬ 
sociated with the dilation of space, the primary reason why 
memory of the initial conditions is lost within a doubling 
time of the scale factor. 

A final (and we believe widely held) misconception is 
that the slope owes to adiabatic processes. This misconcep¬ 
tion arises because the asymptotic 7 —1 is close to the expec¬ 
tation of purely adiabatic evolution, 7 — 1 = 2/3. However, 
intergalactic gas is constantly ascending to a new adiabat 
because of photoheating rather than staying on a single adi¬ 
abat. This ascension is easily noted in the photoheating- 
only case where To^im oc This temporal scaling is 

much different than the scaling a~^ that would occur if 
the evolution were purely adiabatic. After H l reioniza¬ 
tion but before He ll reionization, the heating rate per 


® The ITheuns et al.l lll998l ) approach is equivalent to taking 
the photoheating-only temperature equation (eqn. O and drop¬ 
ping the adiabatic expansion/contraction term owing to inhomo¬ 
geneities (the 2"^^ term on the RHS of eqn. 0 as well as dT/dt 
term, such that the resulting formula equates the photoheating 
rate with the 2HT adiabatic cooling rate. To investigate the 
ITheuns et al.l lll998ll approach in detail, let us drop just the adia¬ 
batic expansion/contraction term owing to inhomogeneities. The 
resulting equation is 

-2/3 _ 2eHirnHia^/^ 

3- “ T77’ 

SkBrihHofdm 

where bars denote the quantity evaluated at A = 1. Solving 
this equation and taking the limit fj ^ fji yields a much flatter 
7 — 1 = 0.3 for collapse histories that have the same n. Figure [l] 
corroborates this result (see th e dA/dt = 0 curve). E quation II27II 
encompasses the one solved in ITheuns et al.l <ll998l l , being more 
general because we have not imposed dT/dt = 0. The approxi¬ 
mation of dropping the d A/dt terms is also made in the He II 
reionization simulations of lCompostella et al.l 1 I 2 OI 3 I I. which may 
explain why their simul ations resulted in some what smaller 7 — 1 
than the simulations of lMcQuinn et ^ (l2009h . 


baryon is 4000(2 -|- abk )~'^^3 A° ® K Gyr“^, and it becomes 
12, 000(2 + ahk)~^Z 3 A° ® K Gyr“^ after He ll reionization 
(Table [ 1 ]). 


7 CONCLUSIONS 


Many studies have noted that the IGM thermal state asymp¬ 
totes quickly to a tight power-law T — A relation. Here we 
showed that this occurs because the photoheating rate is 
nearly constant per driving all cosmic gas to nearly a 

single adiabat in less than a doubling time of the scale fac¬ 
tor. Both cosmic expansion and Compton cooling off of the 
CMB act to accelerate this convergence. That a power-law 
T — A relation arises is a coincidence of how the recom¬ 
bination rate scales with temperature. We showed that the 
final temperature of a gas parcel has almost no sensitivity to 
its previous density evolution, explaining why this relation 
holds so tightly over two decades in density. Furthermore, 
we showed that all relevant cooling processes either reinforce 
the T — A relation established by photoheating or are just 
weak enough not to break the relation at relevant densities. 

We were able to derive the full solution to the temper¬ 
ature in the applicable limit of only cosmic expansion and 
Compton cooling: 


T = 


a? A 


a^Ai 




1/1.7 


where has an extremely weak dependence on the prior 

density evolution of a gas parcel. Further, we showed that 
other cooling processes cause this solution to deviate at < 
10% at A < 10, and so this expression can often be used in 
lieu of a full numerical calculation. 

Cosmological hydrodynamical simulations that employ 
uniform radiation backgrounds show dispersion in the T — A 
relation in excess of simple models based on the Zeldovich 
approximation. This dispersion owes to shocking (and not 
to the collapse history of gas clouds). Such simulations do 
not capture the fluctuations in the gas temperature that 
are an inevitable byproduct of reionization processes. We 
showed that this additional source of fluctuations is erased 
over the same timescale that it takes the mean temperature 
to approach its asymptote. 
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